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Abstract 

Background: Disturbances in the myocardial extracellular volume fraction (ECV), such as diffuse or focal myocardial 
fibrosis or edema, are hallmarks of heart disease. Diffuse ECV changes are difficult to assess or quantify with 
cardiovascular magnetic resonance (CMR) using conventional late gadolinium enhancement (LGE), or pre- or post- 
contrast T1 -mapping alone. ECV measurement circumvents factors that confound T1 -weighted images or T1-maps, 
and has been shown to correlate well with diffuse myocardial fibrosis. The goal of this study was to develop and 
evaluate an automated method for producing a pixel-wise map of ECV that would be adequately robust for clinical 
work flow. 

Methods: ECV maps were automatically generated from T1-maps acquired pre- and post-contrast calibrated by 
blood hematocrit. The algorithm incorporates correction of respiratory motion that occurs due to insufficient 
breath-holding and due to misregistration between breath-holds, as well as automated identification of the blood 
pool. Images were visually scored on a 5-point scale from non-diagnostic (1) to excellent (5). 

Results: The quality score of ECV maps was 4.23 ± 0.83 (m ± SD), scored for n = 600 maps from 338 patients with 
83% either excellent or good. Co-registration of the pre-and post-contrast images improved the image quality for 
ECV maps in 81% of the cases. ECV of normal myocardium was 25.4 ±2.5% (m±SD) using motion correction and 
co-registration values and was 31.5 ±8.7% without motion correction and co-registration. 

Conclusions: Fully automated motion correction and co-registration of breath-holds significantly improve the 
quality of ECV maps, thus making the generation of ECV-maps feasible for clinical work flow. 
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Background 

Cardiovascular magnetic resonance (CMR) using late 
gadolinium enhancement (LGE) provides excellent de- 
piction of myocardial infarction (MI) and focal scar, and 
has become an accepted standard for assessing viability 
[1]. LGE is also useful to detect and characterize fibrosis 
that is "patchy" in appearance associated with non- 
ischemic cardiomyopathies such hypertrophic cardiomy- 
opathy (HCM) [2], Diffuse myocardial fibrosis is more 
difficult to distinguish using LGE since the myocardial 
signal intensity may be nearly isointense and may be glo- 
bally "nulled" thus appearing to be normal tissue [3]. 
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Quantitative measurement of the myocardial tissue lon- 
gitudinal relaxation time constant (Tl) following admin- 
istration of an extracellular Gadolinium-based contrast 
agent is sensitive to increased extracellular volume asso- 
ciated with diffuse myocardial fibrosis but has limitations 
due to a variety of confounding factors [4,5] such as 
gadolinium clearance rate, time of measurement, injected 
dose, body composition, and hematocrit. These factors 
cause a significant variation in post-contrast Tl making 
it difficult to distinguish diseased and normal tissue 
based on absolute Tl values. Pre-contrast Tl varies with 
water content and may be elevated in cases of diffuse 
myocardial fibrosis. Pre-contrast Tl also varies signifi- 
cantly with field strength [6]. Direct measurement of 
extracellular volume (ECV) was initially developed for 
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quantifying the myocardial extracellular fractional distri- 
bution volume [7] and has been proposed as a means for 
detection and quantification of diffuse myocardial fibro- 
sis [3,8-13]. This approach is based on the change in Tl 
following administration of an extracellular contrast 
agent and circumvents the limitation of a single post- 
contrast Tl measurement in detecting a global change in 
Tl. The myocardial ECV is measured as the percent of 
tissue comprised of extracellular space, which is a physio- 
logically intuitive unit of measurement and is independ- 
ent of field strength [6]. ECV has been shown to 
correlate with collagen volume fraction [9,10]. This topic 
is of current interest as a diagnostic tool for non- 
ischemic cardiomyopathies as well as for understanding 
aging processes [9,11-13]. The myocardial ECV may also 
be elevated in edema [14] and may have benefits for 
quantitative assessment as compared to simple Tl and 
T2 measurements. 

Current methods of generating ECV maps require 
cumbersome off-line post processing [13]. The goal of 
this study was to develop and evaluate an automated 
method for producing a pixel-wise map of ECV that 
would be adequately robust for clinical work flow. The 
benefits of pixel-wise mapping over more commonly 
used ROI based measurements is the ability to directly 
visualize the spatial extent and severity of inhomogene- 
ities in ECV in relation to anatomical boundaries. A 
major challenge is the correction of motion that occurs 
due to inconsistent breath-holding during the acquisi- 
tion of images and due to patient movement between 
breath-holds. The proposed method uses a new method 
for motion correction of a series of inversion recovery 
images with widely varying contrasts that is used for Tl- 
mapping [15], and incorporates a new scheme for co- 
registration of pre-and post-contrast Tl-maps. 

Methods 

ECV measurement 

The ECV in the myocardium may be estimated from the 
concentration of extracellular contrast agent in the myo- 
cardium relative to the blood in a dynamic steady state 
[7,9,13]. 



ECV = (1 - hematocrit) Y my ° m myopr \ (1) 
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since the change in relaxation rate AR1 (where Rl = 1/ 
Tl) between pre and post contrast is directly propor- 
tional to the Gd-DTPA concentration, AR1 = y [Gd- 
DTPA] (y = 4.5 L mmol" 1 sec" 1 ). A dynamic steady state 
exists for tissues which have a contrast exchange rate 
with the blood which is faster than the net clearance of 
contrast from the blood [7]. A dynamic steady state 



between the plasma and interstitium may be achieved by 
slow intravenous infusion [3,9], or by imaging 15 min 
following an intravenous bolus administration [12,13] for 
normally perfused myocardium, although 15 min may 
not be adequate for recently infarcted myocardium [16]. 
The bolus method was used in this study since it fits well 
with clinical workflow and permits conventional late en- 
hancement imaging at the desired dose. The ECV for- 
mula (Eq [1]) implies that our myocardial ECV 
measurements include both the intra- and extravascular 
space, and is related to the estimate of extracellular 
extravascular volume fraction (Ve) [8,12] which includes 
additional constant factors. The factor (1 -hematocrit), 
which varies between individuals, represents the blood 
volume of distribution (blood ECV) and converts the 
equation from a partition coefficient calculation to a 
myocardial ECV. 

Image acquisition 

Tl -mapping is based on inverting the magnetization and 
acquiring images at different times during the 
magnetization recovery. The acquisition is ECG triggered 
and all images are acquired at the same cardiac phase in 
late diastole using a Modified Look Locker Inversion Re- 
covery (MOLLI) approach [17]. This method has been 
validated [17,18]. Methods for Tl -measurement such as 
cine inversion recovery (IR) [8,19] which acquire the 
measurements of magnetization at multiple cardiac 
phases may be used for measuring the Tl in the myocar- 
dium but are not suitable for mapping with pixel reso- 
lution. Following respiratory motion correction 
(described below), Tl-fitting of the data at each pixel is 
performed by a non-linear least square fit to the expo- 
nential recovery curve abs(A - B-exp(-TI/Tl*)), where 
Tl is the measured inversion time for each acquired 
image, the absolute value (abs) is used since the images 
are magnitude detected. Tl* is the effective time con- 
stant which includes the effect of the image readout [17] 
related to the desired Tl as T1 = T1* (B/A - 1). This 
"Look-Locker" correction used in MOLLI [17] was 
derived analytically [20] for the case of FLASH readout 
with continuous RF, whereas in the case of MOLLI is 
being applied to gated SSFP readout. This leads to bias 
errors in Tl [21] which are reduced by using a relatively 
low excitation flip angle. Imaging was performed on 
1.5 T Siemens Avanto and Espree scanners (Siemens 
Medical Solutions, Erlangen, Germany), equipped with 
45 mT/m and 200 mT/m/s, and 33 mT/m and 170 mT/ 
m/s gradient systems, respectively. The study protocol 
typically included acquiring Tl-maps for 2 slices (mid- 
ventricular short axis and four chamber long axis views) 
both pre-contrast and approx. 15-20 min following 
intravenous administration of 0.15 mmol/kg Gd-DTPA. 
Both pre- and post-contrast Tl-maps were acquired with 
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the same imaging parameters. The original published 
MOLLI protocol [17] acquired 11 images over 17 heart- 
beats. The MOLLI imaging protocol used in this study 
acquired data at 8 inversion times over an 11 heart beat 
breath-hold at end-expiration with 2 inversions. The ini- 
tial protocol used acquired 3 images acquired after the 
first inversion, 3 heart beat pause, and 5 images acquired 
after the second inversion. The protocol was later modi- 
fied to reduce heart rate variability by acquiring 5 images 
after the first inversion, followed by a 3 heartbeat pause 
and then acquire 3 images after the second inversion. 
Both protocols (3-5 and 5-3) acquired 8 images in 11 
heartbeats. 

Typical imaging parameters were: non-selective inver- 
sion pulse, steady state free precession single shot read 
out with 35° excitation flip angle, field of view 360 x 
270 mm 2 , slice thickness 6 mm, minimum inversion time 
110 ms, inversion time increment 80 ms, matrix 
192x130, voxel size 2.1 x 1.9 x 6.0 mm 3 , TR/TE 2.4/ 

1.0 ms, 7/8 partial Fourier plus parallel imaging factor 2 
with 140 ms readout imaging duration. The maximum 
inversion time was approximately 5 sec at 60 bpm. The 
protocol has subsequently been revised for increased 
spatial resolution using a matrix of 256x144, TR/TE 2.7/ 

1.1 ms (200 ms readout imaging duration), with voxel 
size 1.9 x 1.4 x 6.0 mm 3 for heart rates < 90 bpm, and 
192x130 for hearts above 90 bpm. 

T1 -Mapping and motion correction (MOCO) 

Although patients are instructed to hold their breath 
(typically 11 s) during the MOLLI acquisition of images 
used for Tl -mapping, residual respiratory motion fre- 
quently introduces artifacts due to diaphragmatic drift or 
difficulty with the breath-hold. If uncorrected, this mo- 
tion may lead to errors in the pixel-wise estimation of 
Tl thereby degrading the image quality of the Tl-maps. 
The unique challenge inherent to registration of inver- 
sion recovery images acquired at different inversion 
times is the large variation in image contrast (Figure 1). 



Furthermore, different tissue species with differing values 
of Tl such as blood, fat, normal myocardium, edematous 
myocardium, and infarcted myocardium will experience 
signal nulling at different inversion times. Partial volume 
cancellation at the boundaries between pixels of different 
Tls will occur when the 2 species are out-of-phase 
which further confounds the image registration problem. 
The difference in contrast and appearance between 
images is so large that simple pairwise image registration 
using either intensity or mutual information based 
metrics is inadequate. 

In this study we used the novel motion correction al- 
gorithm described by Xue, et al. [15] based on using 
motion-free synthetic images presenting similar contrast 
to the original acquired inversion recovery images. Add- 
itionally, ECV maps were generated from motion cor- 
rected Tl-maps after co-registration and blood pool 
measurements as described in following sections. Using 
the synthetic image registration approach, the difficulty 
of handling large variations of image contrast through- 
out inversion recovery is bypassed. Motion correction is 
done with an iterative approach based on initial Tl esti- 
mates and the generation of a series of motion free syn- 
thetic inversion recovery images which are used at each 
inversion time for registration with the measured 
MOLLI images. The initial Tl-map is estimated using 
images acquired at the shortest and longest values of Tl, 
which results in a crude estimate of Tl sufficient to 
initialize the iterative process. It was determined that al- 
though the contrast between the short and long Tl 
images is not identical, these images were suitable for 
registration since they avoided the problem of tissue nul- 
ling and partial volume cancellation. From the initial es- 
timate of Tl at each pixel, synthetic images are 
calculated exploiting the known exponential form of in- 
version recovery. Synthetic images are calculated for 
each of the measured images using the estimated values 
of Tl. The synthetic images are all at a fixed respiratory 
position and, importantly, each synthetic image has 
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Figure 1 Pre- and post-contrast image series acquired using a Modified Look Locker Inversion Recovery (MOLLI) sequence used for T1- 
mapping. 
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similar contrast to the corresponding measured image, 
making them suitable for intensity based image registra- 
tion. Due to the non-rigid nature of tissue deformation 
caused by respiratory motion, a non-rigid image registra- 
tion algorithm based on the classical optical flow method 
was to estimate the motion deformation [15]. Image warp- 
ing for both MOCO and co-registration utilized a 
5 th -order B-spline interpolator in order to minimize loss 
of spatial resolution. 

Following motion correction, the Tl map is generated 
by a pixel-wise curve fitting using a three parameter sig- 
nal model for MOLLI inversion recovery. The downhill 
simplex minimization algorithm was used [15]. At each 
pixel the sum-of-squares of the residual errors of the 
measured data to the exponential model fit was calcu- 
lated as a goodness of fit map. 

Co-registration 

Pre- and post-contrast image series are acquired in separ- 
ate breath-holds which are typically 15-30 minutes apart. 
Even small differences in respiratory position or changes 
in patient position due to movement will cause significant 
misregistration of the images. Co-registration of the pre- 
and post-contrast images is performed (Figure 2) using a 



non-rigid image registration applied to the longest inver- 
sion time image of the already motion corrected images. 
In this way, the magnetization of the long inversion time 
images is almost completely recovered and the image 
contrast between pre- and post-contrast series will be 
very similar thereby facilitating intensity based image 
registration methods. Co-registration is achieved by ap- 
plying the estimated motion deformation field to the 
post-contrast series. Rather than simply applying the de- 
formation (warping) to the Tl-map or to the already mo- 
tion corrected image series, a composite deformation 
using the original respiratory motion fields of the post- 
contrast series combined with the deformation from the 
co-registration step is applied a single time to avoid an 
unnecessary loss of spatial resolution. The final Tl-map 
for the post-contrast series is calculated by pixel-wise 
fitting. 

The complete processing workflow was fully automated. 
Once the user specifies the study, the software retrieves all 
images for appropriate paired series of pre and post contrast 
MOLLI acquisitions at the same slice location, performs all 
processing, and produces Tl-maps and ECV maps in dicom 
format which may be pushed to the PACS. The processing 
time in the current implementation is approximately 45 sec 
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Figure 2 Schematic flow chart for automatic generation of ECV maps by respiratory motion correction of pre- and post-contrast 
inversion recovery image series, co-registration of each series with each other, and blood pool segmentation. Note that color scales are 
different for pre and post contrast examples. 
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per ECV series including pre- and post-contrast Tl- 
mapping with MOCO and co-registration, and blood 
pool segmentation. 

Blood pool segmentation 

Calculation of the ECV (Eq [1]) requires knowledge of 
the change in blood relaxation rate (ARl blood ). An auto- 
matic method determining a blood region was developed 
based on the pre-contrast Tl-map. The pre-contrast Tl 
of blood is in the range 1300-1800 ms which is greater 
than the Tl of most tissues. A binary mask was created 
from the pre-contrast Tl map using a threshold of 
1250 ms. Isolated pixels which might be attributed to 
noise were removed, and the resultant mask was 
"eroded" to remove edge pixels that might correspond to 
blood-tissue boundaries with intermediate values due to 
partial volume effects (Figure 3). The values for pre- and 
post-contrast blood Tl were calculated as the median of 
all values of the Tl-map identified by the mask. The 
post-contrast Tl-map was, as a result of the co-registra- 
tion, located in the same position as the pre-contrast Tl- 
map. The change in blood relaxation rate was subse- 
quently calculated as ARl b i ood = l/Tl post - 1/Tl pre . 
Automatic calculation of blood Tls were compared with 
values from manually drawn ROIs in N = 50 subjects. 

Scoring 

ECV maps were scored for overall quality based on a 5 
point scale. N = 600 maps corresponding to 338 patients 
imaged consecutively were scored. The rating was based 
on image artifacts, delineation of myocardial blood pool 
borders, and homogeneity of myocardial signal. A score 
of 5 (excellent) was assigned when the myocardial bor- 
ders were crisp and distinct from the blood pool with less 
than or equal to 1 pixel of partial volume border and no 
noticeable geometric distortion. The myocardial signal 



had to be homogeneous (uniform ECV) in areas assessed 
to be normal tissue (i.e., no focal scar on late enhance- 
ment). A score of 4 (good) was assigned when one of the 
above criteria was absent, i.e., there were 1-2 pixels of 
representative myocardial ECV with slightly increased 
border between blood and myocardium, or geometric dis- 
tortion were present but not impacting the ECV values, 
or there were small regions (10%) of inhomogeneity in 
the myocardium. A score of 3 (fair) meant that that 2 cri- 
teria under the excellent category were absent, i.e., only a 
thin stripe of myocardium was present from which to 
sample the representative ECV, or geometric distortion 
was present, or multiple regions had inhomogeneity 
deemed to be artifactual but affected less than 50% of the 
myocardium. A score of 2 (poor) meant that inhomo- 
geneities or geometric distortion affected greater than 
50% of the myocardium making it difficult to distinguish 
regions of normal myocardium from noise or artifacts, 
but that ECV sampling was still possible in isolated areas. 
A score on 1 (non-diagnostic) meant that ECV could not 
be measured reliably anywhere in the myocardium. 
Examples of ECV maps in each scoring categories are 
shown in Figure 4(a). It is important to note that in for 
certain pathologies, the partial volume at the blood myo- 
cardium border may be mistaken for pathology or vice- 
versa, as emphasized in the Discussion. For the purpose 
of the scoring the quality of maps, the readout involved a 
judgment based on comparison of all the maps including 
with and without motion correction and/or co-registra- 
tion, as well as LGE images and patient history. 

Performance of the motion correction and co- 
registration of Tl-maps were scored separately in addition 
to the overall ECV map quality. These sub-scores were 
assessed as follows. Pre-contrast Tl-maps were compared 
before and after MOCO, and were assessed as to whether 
the MOCO Tl-map was better, the same, or worse than 
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Figure 4 Quality of ECV maps was scored on a scale of 1-5 for N = 600 maps across 338 subjects. ECV map images are provided as 
examples for each quality category. 



the Tl-map without MOCO. Post-contrast Tl-maps with 
and without MOCO were similarly scored. The basis of as- 
sessment was the improvement of myocardial borders and 
reduction of artifacts. Performance of the co-registration of 
pre- and post-contrast data was similarly scored according 
to better, same, or worse by comparing the ECV maps but 
with and without the step of co-registration. 

Human subject protocol 

ECV imaging was performed between Sept 2009 and 
May 2011, 952 patients referred for CMR assessment of 
known or suspected heart disease prospectively under- 
went ECV imaging. Hematocrit was measured from a 
venous blood sample drawn just prior to the CMR study. 
This study was approved by the local Institutional Review 
Boards of the National Heart, Lung, and Blood Institute 
and Suburban Hospital, and all subjects gave written 
informed consent to participate. A total of 1680 ECV 
maps were generated by automated off-line processing as 
described above (728 studies imaged matching short and 
long axis views both pre- and post-contrast, and the 224 
remaining studies only acquired a single slice). Scoring of 
ECV map quality was performed on 600 maps (338 sub- 
jects) acquired consecutively, and statistics for automatic 
blood measurements were made based on the complete 
set. Additionally, cine MRI of cardiac function and phase 
sensitive inversion recovery (PSIR) late Gadolinium en- 
hancement imaging [22] was performed on all subjects. 



ECV measurement 

For the purpose of establishing the value of ECV in 
"normal" myocardium, subjects were categorized as 
"normal" if they were without known significant sys- 
temic illness, demonstrated normal biventricular systolic 
function, chamber size, and wall thickness, absent sig- 
nificant valvular dysfunction, were assessed clinically as 
LGE negative and did not have apparent edema as evi- 
denced by elevated pre-contrast Tl [23]. Subjects with 
diabetes or hypertension were excluded. In order to es- 
tablish ECV values for normal myocardium, ROIs for 
normal subjects were drawn in a mid- wall region of the 
myocardium to avoid contamination from partial vol- 
ume with the blood pool. Values were measured on 
ECV maps with and without motion correction/co- 
registration. 

Statistics 

Data are presented as mean and standard deviation. Dif- 
ferences between groups were tested by the paired or 
unpaired £-test as appropriate. Differences in variance 
were tested by the F test. Statistical significance was 
defined as p < 0.05. 

Results 

Blood measurements 

The measured hematocrit was 41.2 ±4.8% (N = 952 
subjects) and ranged from a minimum of 18% to a 
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Figure 5 Performance of motion correction of Tl-maps and co-registration for ECV-maps was scored as to whether or not there was 
visually detectable improvement or degradation. 



maximum of 65%. A comparison was performed be- 
tween the estimates of blood pool Tl made by the 
automatic blood pool segmentation and manually 
drawn ROIs (Figure 3) for N = 50 subjects. The correl- 
ation between Tl -estimates using manual vs auto- 
matic blood pool segmentation was r = 0.95 (P < 0.001) 
for pre-contrast and r = 0.99 (P < 0.001) for post- 
contrast acquisitions. The pre-contrast Tl was 
1473.2 ±61.8 ms and 1478.1 ±64.7 ms (p = 0.07), using 
automatic and manual measurement respectively. The 
post-contrast Tl was 321.9 ±43.5 ms and 



320.9 ±42.6 ms (p = 0.10), using automatic and manual 
measurement respectively. 

Quality scores 

The quality of ECV maps was 4.23 ±0.83 based on 
N = 600 maps corresponding to 338 unique, consecutive 
subjects. This corresponded to 271 excellent, 225 good, 
87 fair, 12 poor, and 5 non-diagnostic (Figure 4(b)). Tl 
maps with motion correction were assessed to be better 
(14%), the same (85%), and worse (<1%) for pre-contrast 
image series as compared to without correction, and 
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Figure 6 Example case of post-contrast Tl-maps acquired with significant respiratory motion which is significantly improved by 
MOCO, particularly in the septal region (arrows). Residual errors from exponential model fit (arbitrary signal intensity units) are plotted versus 
inversion time for a pixel at the edge of the septum (left) and a pixel at the center of the septum (right) before (non-MOCO) and after MOCO. 
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Figure 7 Example case for which co-registration of data between pre- and post-contrast Tl-maps acquisitions was critical. In these 
examples, MOCO co-registration of the respective maps markedly improved the ECV map, eliminating an artifactually high ECV in the inferior 
septum (arrows) in the non-co-registered maps. MOCO did not improve the quality of individual (pre- and post-contrast) T1-maps for this 
example, however MOCO + co-registration was significantly improved. 



similarly were better (17%), the same (82%), and worse 
(<1%) for post-contrast (Figure 5). A small number of 
instances arose in which the MOCO resulted in signifi- 
cant geometric distortion leading to mis-registration 
which worsened the overall score. An example of a case 
where MOCO improved the Tl-map is shown in Figure 6 
comparing the non-corrected (left) and corrected (right) 
map, and corresponding maps of Tl -fitting error. The 
Tl -fitting error is plotted versus the inversion time for a 
representative edge pixel and a midmural pixel to further 
illustrate the improvement after motion correction. Co- 
registration of the pre-and post-contrast images 
improved the ECV maps in 81% of the cases (example 
shown in Figure 7), were the same in 17%, and were 
worse in 2% of the cases. 

ECV measurements 

Patient characteristics and risk factors for the 62 subjects 
used to define the range of normal values are listed in 
Table 1. ECV of normal myocardium was 25.4 ±2.5% 
using motion correction and co-registration measured in 
the mid-wall for 62 subjects which were LGE negative 



Table 1 Patient characteristics and risk factors for the 62 
subjects used to define normal values 



Age 


43.6 ±17.4 


Gender 


30 male, 32 female 


Family History of CAD 


10 


Smoking 


1 current / 13 ex-smokers 


High Cholesterol 


14 


Hypertension 


0 


Diabetes 


0 


BMI (kg/m 2 ) 


26.5 ±4.6 


Creatinine (mg/dL) 


0.88 ±0.1 9 


LVEF (%) 


61. 7 ±5.9 


LVEDV (mL) 


143.9 ±32.9 


LVEDV index (mL/m 2 ) 


76.8 ±12.2 


LV mass (g) 


87.4 ±26.2 


LV mass index (g/m 2 ) 


46 .2 ±9.5 


Wall Thickness (mm) 


7.5 ± 1 .8 



fiM/=Body Mass Index, LV-\eft ventricular, £F=ejection fraction, EDV-endi- 
diastolic volume. Wall thickness was measured in the basal anteroseptum. 
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Figure 8 Measured ECV values (n = 62 subjects) for normal 
myocardium measured in the mid-wall of the myocardium with 
and without motion correction and co-registration showing 
that motion correction significantly reduces variability 
(p< 0.001). Box and whisker plots show median, 25 and 75 
percentiles, and range. 



and appeared clinically normal Without motion correc- 
tion or co-registration the apparent ECV values were 
30.8 ±7.8% (P< 0.001 compared to motion correction 
and co-registered). The ECV measurements without mo- 
tion correction and co-registration had a greater vari- 
ation due to contamination with blood or other tissue 
(7.8% vs 2.5%, P < 0.001, Figure 8). The ECV of normal 
myocardium with motion correction had a normal range 
(mean ± 2SD) which was 20.4-30.4%. 

Discussion 

Benefits of pixel-wise ECV mapping 

ECV mapping appears promising to complement LGE 
imaging in cases of more homogenously diffuse myocar- 
dial disease states which affect the myocardial extracellu- 
lar space. The ability to display ECV maps in 
quantitative units that may be interpreted on an absolute 
scale offers the potential for simplified detection and 
measurement of the extent of abnormalities affecting the 
myocardial ECV. There is a continuous spectrum of 
spatial heterogeneity between diffuse and focal fibrosis. 
The use of pixel-wise mapping offers a new tool to better 
assess the heterogeneity of the myocardial tissue and to 
measure ECV of the diseased regions. 

Limitations 

The ability to distinguish partial volume border pixels and 
correctly determine whether pixels are contaminated by 
partial volume effects, or are true pathophysiology such as 
subendocardial or subepicardial fibrosis, is a residual and 
important issue. This issue is not only important in terms 
of the visual readout, but may introduce biases into quan- 
titative measurements which become particularly signifi- 
cant in the context of more subtle diseases. One approach 
to reducing the bias in measurements is to restrict the 
measurement to the mid-wall region, although one must 



exercise caution for subjects with thin walled myocardium 
or patients with thin rims of subendocardial fibrosis. 

The proposed approach to generating ECV maps 
resulted in excellent to good image quality maps in 83% of 
the cases. A number of factors contributed to cases of 
poor quality images. In some cases patient movement 
resulted in the slices being substantially different in ap- 
pearance despite the fact that the same slice position was 
prescribed. Co-registration is by design only capable of 
correcting for in-plane motion and cannot compensate for 
through plane motion, although it may appear to partially 
compensate. In other cases, it is clear that the pre- and 
post-contrast images are at substantially different cardiac 
phases due to significant change in heart rate, despite 
adjusting the trigger window. Both of these factors might 
be improved at the time of acquisition by more careful in- 
spection but this would complicate the workflow. 

The interpretation probably can be improved by exam- 
ining the pre-and post-contrast Tl-maps for consistent 
appearance of position and cardiac phase as a quality 
control step. Furthermore, it is expected that the use of 
ECV maps will be in conjunction with pre-contrast Tl- 
maps and conventional late enhancement images, and 
not read independently. As such, the ability to assess and 
quantify diffuse fibrosis in ECV images should provide 
increased diagnostic confidence since LGE image quality 
is generally excellent and will not share the same arti- 
facts. Finally, the use of residual error maps (Figure 7) 
offers a potential as another quality metric to increase 
the confidence for a given region. In other words, if the 
error map indicates that the motion correction is high 
quality (low error), then the ECV values are more 
accurate. 

Partial volume effects affect the measurement of ECV. 
Partial volume effects are a result of limits in in-plane 
spatial resolution as well as slice thickness and are 
observed at tissue borders such as for pixels between 
blood and myocardial tissue. Partial volume effects also 
affect the measurement of fine structures such as MI 
contributing to a larger variability. These same factors 
are also a limitation in LGE. These effects are less of a 
factor for measurements of globally diffuse disease. 

The issue of achieving a "dynamic steady state" using a 
bolus [12,13] rather than slow infusion [3,9] has been 
studied by several researchers, however, the accuracy of 
this assumption and the corresponding impact on quan- 
titative ECV measurements is unknown. In particular, it 
has been noted that steady state may take longer to 
achieve in acute myocardial infarction [16], therefore fur- 
ther work is required to investigate this limitation. 

Despite these issues, ECV maps may be readily incor- 
porated into the clinical workflow and may provide diag- 
nostic information where other methods are limited. It is 
not intended to replace existing methods such as LGE 
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which are excellent at depicting focal lesions, but rather 
to be used in concert with other techniques. When ECV 
is abnormally elevated, it may not be clear whether this is 
due to fibrosis or edema which may be either diffuse or 
focal In such instances, pre-contrast Tl or T2 maps, in 
addition to patient history and contextual imaging clues 
like signs of heart failure, may help to differentiate these 
mechanisms. 

The presence of intra- myocardial fat may affect the 
estimates for ECV and has not been studied. In cases of 
fat with low fat fraction, voxels with partial volume of fat 
and myocardial tissue may lead to errors in Tl -mapping 
which assume a mono-exponential model for single Tl- 
species. Furthermore, fat is not in fast chemical exchange 
with the Gadolinium contrast agent and will not experi- 
ence the same Tl shortening. Multi-echo fat water sepa- 
rated inversion recovery may be used for late 
enhancement [24,25] but is challenging to incorporate 
into a Tl -mapping sequences. Fat water separated im- 
aging may also be used in the combined assessment to 
identify the presence of intra-myocardial fat. 

In this study, the generation of ECV maps was per- 
formed off-line and retrospectively. In order to wholly 
automate the generation of ECV images on the scanner 
at the time of imaging, hematocrit would either have to 
be entered manually or automatically queried from a 
clinical database. These options are potentials for im- 
provement which remain to be explored in order to 
make the generation of ECV images wholly automatic on 
the scanner as a part of the clinical work flow. While it 
is possible to compute a partition coefficient map with- 
out using the value of hematocrit [26], the variability in 
hematocrit (18.3% to 65% in this series) makes the parti- 
tion coefficient less useful for measuring subtle changes 
in ECV expected for diffuse fibrosis. 

The method for automatic segmentation of the blood 
pool was effective for mid-ventricular short axis or long 
axis views acquired in this study. For volumetric cover- 
age which would include apical slices with a reduced size 
of the region of blood, it might be necessary to adapt this 
method to use blood values acquired at other slices with 
the proviso that they must be acquired in close time 
proximity such that the blood Tl remains the same. 

Tl -mapping accuracy and precision 

The accuracy and precision of Tl -mapping using the 
MOLLI method is affected by a number of factors includ- 
ing the protocol, the sequence design, patient specific fac- 
tors related to heart rate and blood flow, and system 
adjustments such as center frequency and shim. The 
MOLLI method for Tl -mapping is widely use and a num- 
ber of adaptions to the protocol have been introduced to 
shorten the acquisition time. Further optimization the se- 
quence and of pre- and post-contrast Tl -mapping 



protocols based on the MOLLI strategy is the subject of 
ongoing research. 

A shortened acquisition acquiring 7 images in 9 heart- 
beats (known as ShMOLLI) [27] may be used to further 
reduce the breath-hold duration and achieve a heart rate 
independent Tl measurement by discarding samples for 
tissue with longer Tl. The 11 heart beat protocol used in 
this study achieves a reliable Tl measurement by using a 
longer recovery between inversions without eliminating 
data. At longer Tls the 9 heart beat ShMOLLI protocol 
uses only 5 images for Tl -fitting, whereas the proposed 
motion corrected 11 heart beat protocol uses all 8 
acquired images. The increased number of images 
reduces the fitting error. 

In terms of the sequence design, we have recently 
determined that imperfect inversion efficiency leads to 
an underestimation of Tl, representing the largest 
source of systematic error. This error may be greatly 
reduced by means of an improved inversion pulse design 
as well as by applying an appropriate correction. These 
corrections were not available at the time of this study. 
As a result, the normal values for Tl are underestimated 
by as much as 5-10%. However the systematic bias in 
Tl -measurement results in only a small overestimation 
of ECV, estimated on the order of 1% or less in ECV 
units. This does not alter the main conclusions of this 
present study on the utility of ECV maps. 

Conclusions 

While LGE is an excellent method for detection of focal 
scarring as encountered in myocardial infarction and nu- 
merous other cardiomyopathies, ECV mapping appears 
promising to augment LGE imaging in cases of more 
homogenously diffuse disease. The fully automated im- 
plementation lends itself to incorporation of ECV map- 
ping directly into a clinical workflow. Motion correction 
and co-registration of breath-holds have been demon- 
strated to significantly improve the image quality of ECV 
maps. 
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